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摘 要 : 无 线 传感器 网 络 的 最 大 干扰 最 小 化 问题 可 以 被 描述 为 已 知 平面 上 nn 个 点 的 位 置 及 发 射 半径 的 国 值 ， 要 求 设 
置 发 射 半径 使 得 任意 点 被 其 他 点 发 射 范围 覆盖 的 最 大 数量 达到 最 小 。 为 了 有 效 求解 该 问题 ， 提 出 了 一 种 新 的 贪心 算 
法 一 一 缩 边 算法 。 不 同 于 已 有 算法 的 构图 方式 ， 该 算法 是 通过 采取 缩 边 的 方式 ， 构 造 网 络 通信 拓扑 图 ， 并 结合 了 操 
作 系 统 中 批量 处 理 的 思想 对 贪心 算法 进行 了 加 速 ， 缩 短 了 算法 的 运行 时 间 。 通 过 实验 验证 ， 该 算法 相 比 于 已 有 算法 
在 随机 产生 的 算 例 上 能 产生 更 优 的 解 。 
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Edge shrinking greedy algorithm for minimizing interference in wireless networks 
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(1. School of Computer, South China Normal University, Guangzhou 510631, China; 2. School of Mathematics & Statistics, 
Shaoguan University, Shaoguan Guangdong 512005, China) 


Abstract: The problem of minimizing the maximum interference in wireless sensor networks can be described as follows: 
given n points with their positions on the plane and thresholds of transmission ranges , for each point, the maximum number 
of other points whose transmission ranges cover this point is minimized. To effectively solve this problem, this paper 
proposed a new greedy algorithm called edge shrinking algorithm. The algorithm constructs a desirable network 
communication topology by shrinking edges, which is different from existing algorithms and is accelerated by incorporating 
the idea of batch processing in operating system to shorten its running time. The result shows that this algorithm can 
produce better solutions than existing algorithms. 
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在 平面 上 计算 无 线 传感器 网 络 的 最 大 干扰 最 小 化 并 保持 
网 络 的 连通 已 经 被 证 明 是 NP 难度 的 向 ， 因 此 要 在 多 项 式 时 
| 

多 


0 引言 


FF 线 传 感 器 网 络 (wireless sensor network, WSN) 是 由 大 间 内 找到 精确 的 解 是 困难 的 。 文 献 中 求解 该 问题 的 算法 并 不 
量 的 静止 或 移动 的 传感器 以 自 组 织 和 多 跳 的 方式 构成 的 无 线 ， 典 型 的 有 如 下 几 种 : 2008 年 ，Halldorsson 提出 的 结合 

网 络 ， 以 协作 地 感知 、 采 集 、 处 理 和 传输 网 络 履 盖 地 理 区 域 计算 几何 和 s- 网 络 理论 的 近似 算法 中， 其 最 大 干扰 是 405， 其 
'- 内 被 感知 对 象 的 信息 ， 并 最 终 把 信息 发 送 给 网 络 的 所 有 者 。 中 4 是 发 射 半径 相同 的 网 络 中 节点 最 大 度 。 同 年 又 有 学 者 证 
FL 无 线 传感器 网 络 具 有 的 众多 类 型 的 传感器 ,可 探测 包括 地 震 、 明了 任何 能 够 产生 连通 图 的 传输 半径 函数 都 可 以 修改 为 具有 
电磁 、 温 度 、 湿 度 、 噪 声 、 光 强度 、 压 力 、 土 壤 成 分 、 移 动 以 下 特点 的 半径 函数 ,其 产生 的 拓扑 结构 干扰 基本 不 变 , 且 所 
物体 的 大 小 、 速 度 和 方向 等 周边 环境 中 多 种 多 样 的 现象 (1。 有 传感器 节点 的 发 射 半 径 小 于 Rmin Rmin 为 保持 网 络 强 连通 
无 线 传感器 网 络 一 般 部 署 在 人 们 不 方便 直接 控制 的 地 方 ， 的 最 小 可 能 半径 四 。2009 年 ， 又 有 人 提出 基于 每 个 节点 与 最 
而 且 传感器 节点 一 般 使 用 的 是 电池 ， 很 难 重新 充电 或 者 更 换 ”近邻 居 关 系 的 最 近邻 算法 四。2010，Aslanyan 等 人 提出 了 一 
电池 ， 因 此 能 源 成 为 制约 无 线 网 络 服 务 时 限 的 关键 因素 之 一 种 基于 部 分 成 员 和 集合 履 盖 的 近似 算法 中。2014 年 ， 有 人 基于 
DB]。 节 省 节点 能 源 的 一 种 方法 是 减少 附近 节点 同时 传输 信和 号 Prim 算法 提出 了 一 种 时 间 复 杂 度 为 O(n?) I e Pe S SETA, 


所 产生 的 干扰 。 之 后 又 有 人 尝试 在 最 佳 邻 算 加 入 了 后 悔 机 制 0 轨 ， 在 适当 增 
前 减少 干扰 的 方法 有 频率 分 配 多 路 技术 、 时 间 分 配 多 ”加 时 间 的 前 提 下 ， 一 定 程度 上 提高 了 算法 的 精确 度 。 

路 技术 和 拓扑 控制 三 种 。 前 面 两 种 方法 虽然 都 可 以 有 效 地 减 本 文 考虑 了 用 拓扑 控制 方法 来 减少 平面 中 无 线 传感器 网 
少 干扰 ， 但 是 要 求 较 高 的 技术 和 成 本 ， 而 拓扑 控制 是 通过 给 。” 络 中 最 大 干扰 的 最 小 值 ， 提 出 了 一 种 不 同 于 已 有 算法 构图 方 
每 个 节点 分 配 一 个 合理 的 发 射 半径 并 形成 一 个 连通 的 拓扑 ， 。 式 的 新 的 贪心 算法 ， 并 在 大 量 的 随机 算 例 上 验证 其 有 效 性 。 
并 且 减 小 半径 的 某 个 目标 函数 使 网 络 的 最 大 干扰 最 小 因此， DD 

在 无 线 传感器 网 络 中 ， 拓 扑 控制 成 为 了 减少 干扰 的 重要 方法 1 ”问题 模型 以 及 相关 算 ; 

之 一 内， 也 是 近年 来 的 一 个 研究 热点 。 针对 无 线 通信 中 减少 干扰 的 问题 现 已 提出 多 个 不 同 的 模 
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4o 4, Si 无 线 网 络 干 


。 其 中 有 一 种 是 以 发 送 者 为 中 心 的 模型 
F 扰 。 另 一 种 以 接收 者 为 中 心 的 模型 中 ， 了 


4， 在 此 模型 中 了 


点 ”的 干扰 定义 为 传输 范围 
2 以 发 送 者 为 中 心 的 模型 只 考虑 
在 接受 者 上 ， 模 型 不 能 反映 真 3 


ta 


于 将 无 线 传感器 网 络 建 模 为 图 


FAT 1 
SK HJ 


FP 心 的 模型 ， 通 过 拓扑 控制 


以 接受 者 为 中 心 模型 的 无 线 传感器 网 络 中 干扰 最 小 化 问题 的 


论 的 知识 来 描述 问题 模型 ， 因 此 需要 了 


目 关 文献 49。 下 面 给 昌 
的 定义 并 介绍 最 小 干扰 
的 无 线 传感器 网 络 可 以 


4 络 ” 是 通用 的 。 关 于 图 的 相关 概念 可 
上 相关 的 概念 ， 然 后 给 出 形式 化 
化 模型 。 


合 中 的 每 个 点 i 都 有 


径 函 数 r(i) € [0,RG)]» 


作出 以 下 定义 。 


的 传输 关系 i 


大 点 u,v EV， 经 由 传输 半径 函数 x 


述 为 节点 的 集合 ， 在 全 
住 坐标 (x(i)，y( 让 )， 并 且 满 足 其 
中 RQ) 是 该 节点 的 半径 阔 值 。 因 此 ， 
于 传输 半径 函数 + 所 导出 的 所 有 节点 的 传输 关系 可 以 建 模 
到 G,-«V,Er». XF V PRW u 和 vw， 边 (u,v) 
存在 当 且 仅 当 rM r(v) P T. u 和 vw 的 欧 氏 距离 d(w,v)， 即 
r(v)>d(u,v)。 本 文 为 了 方便 起 见 ， 将 G,=<V,E,> 
叫做 传输 半径 函数 的 导 
定义 1 对 六 中 任意 的 7 


UN 7 


ChinaXiv 合 作 期 刊 


光 最 小 化 问题 的 一 种 缩 边 贪 心算 法 第 37 卷 第 5 期 


它 在 图 G; 中 距离 最 远 的 邻居 的 距离 作为 节点 v 的 传输 半径 ， 
最 终 的 得 到 的 导出 图 和 导出 图 G, 有 相同 的 干扰 。 因此， 可 以 
仅 从 节点 v 与 其 在 导出 图 Ga 中 可 到 达 的 邻居 节点 的 不 同 距 
离 中 来 选择 v 的 传输 半径 。 通 过 这 样 的 选择 方式 ， 节 点 v 的 
和 专 输 半 径 即 导出 图 G; 中 vv 与 可 到 达 的 距离 最 远 的 邻居 的 距离 。 
MIP 已 经 被 证 明 是 NP 难度 的 ， 即 使 对 于 寻找 所 有 节点 
都 在 一 条 直线 上 这 种 特殊 的 MIP 的 一 个 最 优 实例 也 是 很 困 
为 ， 尽 管 难度 还 没有 被 证 明 出 来 。 对 于 二 维 网 络 的 特殊 情 
如 网 格 网 络 ， 这 个 问题 的 难度 也 是 NP 难度 的 。MIP 难 
源 于 其 组 合 结构 ， 以 至 于 从 许多 直观 的 特性 ， 如 最 小 度 、 
最 近 的 邻居 、 最 好 的 邻居 等 特点 , 都 不 能 产生 令 人 满意 的 解 。 
寻 此 可 以 考虑 设计 贪心 算法 来 解决 这 个 问题 。 

1.2 干扰 最 小 化 问题 的 典型 算法 
前 已 有 不 少 学 者 研究 了 干扰 最 小 化 问题 ， 并 提出 了 解 
决 该 问题 的 相关 算法 ， 本 文 主要 介绍 了 介绍 如 下 几 种 算法 。 
1.2.1 最 近邻 算法 (NNA) 
贪心 算法 的 基本 思路 是 从 问题 的 某 一 个 初始 解 
一 步 地 进行 ， 根 据 某 个 优化 测度 ， 每 一 步 都 要 确保 
部 最 优 解 。 在 干扰 最 小 化 问题 中 ， 可 直观 地 认为 节 
半径 越 小 则 网 络 中 的 干扰 越 小 。 最 近邻 算法 正 是 基 
直观 特性 ， 并 且 结 合 了 最 小 生成 树 的 思想 ， 通 过 调 
邻 最 近 的 节点 的 传输 半径 来 连通 当前 网 络 中 不 连通 的 两 个 连 


发 一 步 
获得 局 
的 传输 
这 样 的 


RE H r ORE EE 


5 UKE G=<V,E>, W(u,)E 
-存在 当 且 仪 当 w 和 vw 的 传输 半径 rM r(v) 不 小 于 
的 欧 氏 距离 d(u,v), BB r(u)ed(u,v) B. r(v)ed(u,v). 


F 它们 之 


G,=<VB> 称 为 传输 半径 函数 r+ 的 导出 图 。 
G,-«V, E> PESA v 的 干扰 度 定义 为 
EF 径 能 覆盖 到 节点 v 的 图 中 其 点 的 个 数 ， 即 

1.0v) {u eV \{v} |ru) 2 d(u,v)}| 
| G,-V,E,7 BY LEE 1,58 3C Er Br A 


Lv) Xo 


y 


8 1 中 点 v 的 干扰 是 3, 导出 


I, = max(1,(v)|veV] 


DS 


CK UR AR 


Fig. 1 Example diagram of maximum interference 


一 点 i 都 有 一 个 满足 以 下 条 


基本 的 要 求 ， 可 以 有 下 面 的 问题 模型 
F 扰 化 问题 MP) :给 定 平 十 
标 可 表示 为 CGOD)，?y(D)。 对 二 中 的 任意 
牛 的 传输 半径 阔 值 函数 RC) € [0,1] 


n 
p 


导出 图 Gr=<V,Er> 连 通 
足 导出 图 G,=< 人 E> 是 
从 上 文 的 定义 可 多 


因此 需要 确定 一 个 传输 半径 了 满 


的 且 干 扰 最 小 的 。 


对 任意 一 个 传输 半径 x 都 满足 导出 


WR] 


。 而 且 ， 如 果 重 新 选择 使 节点 v 到 


的 一 个 节点 集 V, 


通 分 量 。 重 复 这 个 步 又， 直至 整个 网 络 连通 。 
1.2.2 最 佳 邻 算法 (BNA) 

最 近邻 算法 在 一 般 网 络 中 的 效果 尚 可 ， 但 是 在 特殊 网 络 
〈 指 数 链 网 络 ) 中 效果 很 差 。 针 对 最 近邻 算法 的 不 足 ， 文 献 


[10] 提 出 了 改进 的 最 近邻 算法 一 一 最 佳 邻 算法 。 最 佳 邻 算 法 
的 主要 思想 是 : 通过 调整 当前 连通 网 络 最 大 干扰 增幅 最 小 的 
节点 对 的 传输 半径 来 连通 两 个 不 连通 的 部 分 。 重 复 上 述 步 又 ， 
直至 整个 网 络 连 通 

最 佳 邻 算法 与 最 近邻 算法 的 时 间 复 杂 度 相当 。 但 是 大 量 
的 实验 结果 表明 ， 最 佳 邻 算法 不 仅 能 消除 最 近邻 算法 在 特殊 


网 络 上 的 次 端 ， 在 普通 网 络 上 的 效果 也 优 于 最 近邻 算法 。 
1.2.3 后 悔 贪 心算 法 (RGA) 

Sunl IÆ 2017 年 提出 了 后 悔 贪心 算法 ， 其 主要 步骤 遵循 
贪心 算法 的 工作 流程 ， 并 在 贪心 算法 的 每 个 迭代 步 尝 试 采用 
后 悔 策略 。 在 贪心 算法 的 每 个 迭代 步 中 ， 如 果 当 前 网 络 的 最 
大 干扰 发 生 了 改变 ， 即 当前 导出 图 的 最 大 干扰 增加 后 ， 将 针 
对 此 时 具有 最 大 干扰 的 节点 用 后 悔 策 略 调整 某 些 节点 的 已 
贪心 算法 确定 的 传输 半径 来 试图 降低 它 的 干扰 ， 使 得 当前 最 
大 干扰 减少 。 
1.2.4 旋转 贪心 算法 (CGA) 
基于 基本 贪心 算法 G 的 旋转 贪心 算法 的 主要 思想 是 在 贪 
心算 法 G 求 得 的 可 行 解 中 选择 一 定 比例 的 元 素 作 为 部 分 解 。 
每 次 将 部 分 解 的 第 一 个 元 素 移 除 , 通过 贪心 算法 G 重新 加 入 
元 素 ， 更 新 部 分 解 。 重 复 这 一 步骤 ， 直 至 迭代 结束 。 最 后 利 
用 贪心 算法 G 得 到 一 个 完整 的 解 。 
1.2.5 基于 部 分 成 员 集合 履 盖 算法 (SCA) 

Kuhn 等 人 在 2005 提出 了 最 小 成 员 集 合 履 盖 问 题 ， 并 将 
其 用 于 解决 蜂窝 网 络 中 的 干扰 问题 。 在 文献 [9] 中 ，Aslanyan 
等 人 在 2010 年 将 最 小 成 员 集合 覆盖 问题 推广 到 最 小 部 分 成 
员 集 合 覆 盖 ， 将 其 用 于 求解 无 线 网 络 干扰 最 小 化 问题 ， 并 从 
理论 上 证 明了 干扰 的 上 限 O(opxiIn?n) RR fo EG O(optxIn?n). 
部 分 成 员 集 合 履 盖 定义 如 下 : 给 定 一 个 n 元 素 集合 SS 
US 以 及 5S 子 集 的 集合 C={C1, C», Ca RP SINSE ó,[Sien, 
Sən, RRE C'c C 且 满 足以 下 两 个 条 件 : 


mm 
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CE Sy 中 的 所 有 元 素 ， 可 表示 为 SS e vec o 
CHE S» 中 的 元 素 被 覆盖 的 最 大 次 数 最 小 ， 即 max M (s C M Iu 
最 小 。 其 中 5; 中 的 元 素 u 被 集合 CUR i EUCROE X 737638 u 4 的 出 度 为 7 e 
的 成 员 值 M(u,C)， 可 表示 为 M=(u,C”)=|CiE Cue Cilo KA) 
然后 将 该 问题 构建 成 一 个 整数 规划 问题 ， 利 用 线性 规划 — 
松弛 技术 解决 ， 因 此 可 以 将 C' 看 做 是 问题 的 一 个 解 ， 对 于 每 BIZ, O22 
一 个 GEC, 4l —^ xE {0,1}, HI x-1e cec. 
Minimize z | 
to xil wes, | 
Ds XSz  ujeS, ! 
x; €(0.1] C;eC $ 
具体 步骤 如 下 : 图 2 缩 边 点 选取 实例 
a) 初 始 化 图 G=<V, Eo», Eo- Q, I-0. Fig.2 Example of selecting shrinking point 
b)/=1t+1， 找 到 问题 中 对 应 的 $1、S。、C， 并 用 线性 规划 图 2 中 直观 展示 了 按 式 (3) 选 取 缩 边 点 的 合理 性 ， 假 设 u 
RC 点 为 当前 最 大 干扰 点 且 最 大 干扰 为 3〈 为 了 方便 理解 ， 将 所 
c) 调 整 C' 中 边 上 的 半径 , EEEn U C', 更 新 图 GF<V,E2>。 有 在 传输 半径 以 内 的 点 都 放 在 同一 直线 上 ) ， 点 4、B、C 分 
中) 重复 步骤 b) 和 c)， 直 至 图 连通 为 止 。 别 为 干扰 u 的 三 个 节点 ， 可 看 出 B 的 出 度 最 大 ， 而 u TC 
EA 的 发 射 范围 的 最 外 围 ， 当 缩小 节点 B 或 节点 C 的 发 射 半径 都 
2 缩 边 贪心 算法 (SEA) 会 使 得 缩 边 的 节点 仅 与 一 个 节点 相连 ， 均 会 破坏 导出 的 连通 
2.1 基本 思想 性 。 根 据 式 (3) 选 取 的 点 为 4， 缩 小 点 4 的 发 射 半径 能 保持 导 
缩 边 贪心 算法 (简称 缩 边 算法 ， 以 下 均 使 用 简称 ) 主要 出 图 的 连通 性 。 Un GRE B ska C 都 会 极 大 地 破坏 导出 
思想 是 从 最 大 传输 半径 RR 的 导出 连通 图 Gr=<V,Er> 开 始 ， 通 图 的 连通 性 。 
过 逐步 缩小 图 中 各 点 的 传输 半径 ， 找 到 干扰 较 小 的 连通 子 图 2.2 算法 描述 
G,=<V,E,>。 算 法 开始 时 通过 将 所 有 节点 的 传输 半径 初始 化 到 在 本 节 中 将 对 缩 边 算法 进行 描述 分 析 ， 并 根据 各 种 算法 
EARB. B ORE], Ev, 即 初始 化 的 导出 图 是 连通 在 相同 实例 上 的 执行 结果 来 比较 不 同 算法 产生 的 解 的 质量 。 
的 ， 在 确保 导出 图 是 连通 的 前 提 下 算法 按照 缩 边 策略 去 减 小 算法 1 给 出 缩 边 算法 的 伪 代 码 。 
各 节点 的 发 射 半径 ， 直 至 所 有 节点 的 传输 半径 都 不 能 再 减 小 算法 1 缩 边 算法 
为 止 。 而 已 有 算法 都 是 通过 初始 化 图 中 各 节点 发 射 半径 为 0 Input :已 知 节点 集 V={1,2, .. .nm 在 平面 上 的 坐标 ， 对 于 每 个 点 1 
开始 ， 然 后 一 步 一 步 地 按照 算法 中 的 策略 扩大 某 些 节点 的 发 EV 的 最 大 传输 半径 阔 值 函数 Ri E [86,1], 满 足 图 Gam (V, Er) 是 连通 
射 半 径 使 得 这 些 点 加 入 到 当前 导出 图 的 连通 分 图 中 ， 即 通过 的 。 
“加 边 ” 的 方式 将 这 些 节点 连通 到 当前 的 连通 图 中 ， 与 之 相 Output: 图 的 传输 半径 向 量 上 = (ra n2, ...，,rn), 使 得 图 Gr=(V,Er) 
比 , 缩 边 算法 又 可 以 看 做 是 不 断 地 对 当前 导出 连通 图 通过 “ 减 连通 且 干 扰 较 小 。 
边 ” 的 方式 以 达到 降低 干扰 的 目的 。 Begin 
于 最 终 所 得 的 导出 图 必须 是 连通 的 ， 所 以 选 点 缩 边 要 1. m=Ri, i€ (1,2,...,n) 
遵循 的 首要 条 件 是 不 能 破坏 导出 图 的 连通 性 。 根 据 这 条 件 选 2. do ue 当前 最 大 干扰 度 点 
择 缩 边 的 节点 ， 直 观 上 所 选 的 点 应 是 出 度 大 且 最 大 干扰 点 应 3. viEV, N(u)—{veV|Irv2d(u,v)} 
在 其 覆盖 圆 的 最 外 围 。 出 度 大 意味 着 该 点 的 连通 性 能 可 能 较 4. While N(u)z0 
好 ， 缩 小 其 发 射 半径 对 导出 图 的 连通 性 影响 不 大 ;而 最 大 干 stn 
扰 点 在 其 覆盖 圆 最 外 围 的 点 则 意味 着 小 幅度 的 缩小 点 的 发 射 20 
半径 即 可 消除 对 最 大 干扰 点 的 干扰 ， 最 小 程度 地 降低 了 对 导 5. if 对 Vv’ 缩 边 后 图 仍 连通 then 
出 图 连通 性 的 影响 。 然 而 “出 度 大 ”和 “最 大 干扰 点 处 在 其 6. do E% rv 
外 围 ” 这 两 点 如 果 仅 仅 单一 考虑 一 个 都 会 出 现 极端 的 情况 ， 7. else 
会 极 大 地 影响 当前 导出 图 的 连通 性 。 为 了 避免 出 现 这 些 极端 8. N(u) -N(u)Mv') 
情况 ， 下 面 综合 考虑 这 两 点 来 选取 缩 边 的 节点 。 为 此 作出 如 9. end if 
下 定义 : 10. end while 
YE u 为 当前 最 大 干扰 点 ， 令 NGCO) 为 对 2 产生 干扰 的 节点 11. while u 的 干扰 减 小 
EE, B Nw)={v EVIr>d(uyv)}, 令 1(u,v) 为 u 在 v 中 按 与 vy 距 12. end while 
离 从 近 到 远 的 排序 号 ， 其 中 vE N(n)， 则 选择 缩 边 的 点 为 end 
vag mx) 3 下 述 定理 给 出 了 缩 边 算法 的 正确 性 及 时 间 复 杂 度 的 结论 。 


定理 1 该 算法 在 O(n3xq) 时 间 内 产生 一 个 传输 半径 函数 

缩 边 点 的 发 射 半径 由 rH NIS du,v)-e， 其 中 为 任意 r HFAB G,=<WE,> 是 连通 的 。 
小 的 正 实数 。 证 明 首先 证 明 算法 的 时 间 复 杂 度 。 该 算法 第 一 步 首先 
要 找到 当前 网 络 中 最 大 干扰 节点 uw， 需要 时 间 O(n); 接着 第 
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二 步 要 找 出 所 有 出 度 中 含有 的 节点 需要 时 间 O(n); 最 后 为 。” 缩 边 算法 的 执行 情况 并 与 已 有 的 算法 作 了 对 比分 析 。 算 例 中 
了 保证 图 的 连通 性 ， 在 每 次 缩小 节点 的 传输 半径 之 后 需要 判 10 个 点 vi,va,...,vio 在 平面 上 的 位 置 如 图 4 所 示 ， 对 应 每 个 点 
导出 子 图 是 否 连 通 需 要 时 间 O(n)。 又 因为 该 算法 需要 的 坐标 分 别 为 (1.11，0.17)(1.64，1.84)(1.92，1.72)(0.46， 
述 直 至 再 也 找 不 到 可 以 缩小 的 节点 为 止 ， 一般 这 个 重 ”0.79)(1.94，0.81)(0.29，1.84)(0.24，1.30)(0.60，1.00)(1.56， 
数 4 会 远 小 于 n， 所 以 该 算法 的 时 间 复 杂 度 为 O(n3xqd)。 0.37)(1.54，2.26)。 

次 证 明 导 出 图 的 连通 性 ， 该 算法 初始 化 图 时 ， 就 确保 

了 图 是 连通 的 ， 并 且 每 一 步 选择 缩小 传输 半径 的 节点 都 必须 > i 
到 的 连通 性 为 首要 条 件 ， 因 此 所 选择 的 点 不 影响 " H 
图 的 连通 性 ， 所 以 该 算法 产生 的 导出 图 G<, E> E XEXÉITT e » 
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am 
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23 算 例 分 析 ! 
f| 1 对 于 上 述 文中 所 提 到 的 一 维 指数 链 网 络 的 实例 ， 
考虑 在 一 维 指数 链 网 络 n=6 的 实例 上 用 缩 边 算法 的 执行 情况 。 
初始 化 图 时 ,所 有 的 节点 都 连接 到 了 与 自己 相距 最 远 的 节点 ， M 
如 图 3(a) 所 示 ， 其 中 每 个 点 的 弧 线 表示 该 节点 的 传输 半径 最 l 
远 能 影响 到 的 节点 。 图 4 节点 在 平面 上 的 位 置 

Fig.4 Position of points on plane 
已 知 上 述 平 面 坐标 的 节点 集 天 142…,10}， 算 例 的 最 大 
T DEFINES "S i in ts E He HERE. MER Ga-(V.ER) AE VERI) RE 
执行 时 进行 初始 化 ， 即 将 每 个 节点 的 传输 半径 设置 为 所 能 连 
(a) 在 6 个 点 的 一 维 指 数 链 网 络 , 缩 边 算法 初始 化 各 点 半 和 接 到 最 远 的 节点 。 初 始 化 示意 图 如 图 5 所 示 ， 图 中 箭头 表示 
产生 的 导出 图 该 节点 所 能 到 达 最 远 的 节点 。 


(a) In one-dimensional exponential chain network with six points, edge 


N 


HX 
Fi 


reduction algorithm initializes derived graphs generated by radius of 


each point (1)-091 
1(2)-0.45 
1(3)-0.91 
1(4)=0.91 

vı vò V3 V4 Vs V6 1(5)-0.89 
1(6)-0.89 
1(7)20.55 
Š r(8)=0.91 
(b) 执行 一 次 算法 迭代 步 后 的 导出 图 r(9)=0.91 
r(10)=0.66 
(b) Derived graphs after performing iteration step of algorithm inc 
W v WV V4 Vs V6 
(c) 在 6 个 点 的 一 维 指 数 链 网 络 ， 缩 边 算法 的 最 终结 果 导 出 区 Q5 各 节点 传输 半径 初始 化 
(c) In one-dimensional exponential chain network with six points, final Fig.5 Initialize transmission ranges of points 


result of edge reduction algorithm is derived 
图 3 缩 边 算法 在 一 维 指 数 链 实例 
Fig. 3 Example of shrinking algorithm in one-dimensional chain 043 
此 时 图 中 6 个 节点 的 干扰 度 都 为 5， 均 为 最 大 干扰 点 。 
最 大 干扰 节点 不 止 一 个 时 ， 按 顺序 先 队 编号 小 的 点 先进 行 缩 v 
边 ， 选 取 的 需要 减少 干扰 的 节点 为 vw， 则 缩 边 节点 


”ton 。 依 此 类 推 ， 接 下 来 需要 减少 干扰 的 节点 


veN(w) 


Vio 


0.54 


分 别 为 va. vss va XIM RJAR TS EXT 973. vos ve 和 v1， 完成 li 

这 些 缩 边 操作 之 后 ， 各 节点 传输 半径 如 图 3(b) 所 示 。 经 过 上 

述 一 系列 缩 边 之 后 ， 节 点 vi 和 节点 ve 都 缩小 了 传输 半径 ,使 0.25 

得 当前 图 的 最 大 干扰 为 4。 

重复 上 述 步 又 ， 直 至 不 存在 可 以 缩 边 的 节点 为 止 ， 最 终 

结果 示意 图 如 图 3(c) 所 示 。 图 3(c) 中 ， vi、v2、v4、vVs 和 ve 的 

传输 半径 都 仅仅 覆盖 到 了 与 自身 相 邻 最 近 的 两 个 节点 《〈 除 端 Yo 

节点 外 ) ， 而 节点 v3 的 传输 半径 则 覆盖 到 了 与 自身 相 邻 最 远 n 

的 一 个 节点 即 ve。 从 示意 图 中 可 直观 地 观察 到 导出 图 是 连通 图 6 最 终 各 点 传输 半径 

的 ， 且 最 大 干扰 度 为 2， 平 均 干 扰 度 为 1.83. Fig.6 Final transmission ranges of points 
例 2 考虑 在 二 维 空间 随机 产生 了 n=10 个 点 的 算 例 上 用 缩 边 算法 每 个 迭代 步 中 节点 传输 半径 变化 如 下 : 


Z 


0.89 


201904.00046v1 
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a)r(1)=0.91—0.89 

b)r(8)-0.910.89 

c)(10)-0.660.43 

d)r(7)-0.550.54 

e)((8)-0.890.25 

£(6)-0.890.54 

最 终结 果 示 意图 见 图 6。 为 了 清晰 起 见 ， 这 里 仅 标 记 
了 经 过 一 系列 缩小 传输 半径 后 各 个 节点 传输 半径 内 所 最 远 
到 达 的 节点 。 从 图 中 各 节点 的 传输 半径 可 知 ， 导 出 图 是 连 
的 ， 最 大 干扰 点 是 有 多 个 分 别 是 vi. va. vs. va. vs. v. vs. 
v9， 其 干扰 值 均 为 2。 
在 该 算 例 中 前 文 提 到 的 “加 边 ” 算 法 在 这 个 算 例 中 都 得 
出 了 相同 的 结果 ， 其 中 最 大 干扰 为 3， 如 图 7 所 示 。 


ER BE EE 


Vio 
0.43 
V6 
v2 
0.54 0457 3 
v7 0.47 0.91 
Vg 
Va Vs 
0.5 

0.89 

V9 


图 7 加 边 ” 算 法 最 终结 
Fig.7 Result of adding edge algorithm 

将 图 6 与 7 对 比 可 看 出 ， 图 6 中 唯一 的 区 别 在 于 缩 边 算 
法 中 设置 x(8)=0.25， 而 “加 边 ” 算 法 中 设置 x(8)=0.47， 从 而 
导致 最 终 最 大 干扰 值 不 同 。 因 此 在 这 个 实例 上 缩 边 算法 所 得 
的 解 比 “加 边 ” 算 法 得 到 的 解 更 优 。 
2.4 实验 结果 

实验 环境 如 表 1 所 示 。 
AY 实验 环境 说 明 


Table 1 Experimental environment description 
Type Specific information 
OS Microsoft Windows 7 Pro version 64bit 
CPU Intel(R) Core(TM) i5-2410M CPU @ 2.30GHz 
Memory 6GB 
Disk HHD 512GB 
IDE Codeblocks 13.12 


ChinaXiv 合 作 期 刊 
第 37 卷 第 5 期 


表 2 各 个 算法 所 得 解 的 平均 目标 值 
Table 2 Average result of algorithms 


(n,d) NNA BNA RGA CGA SCA SEA 


(100,6) 5.01 4.62 4.40 4.41 4.74 4.02 
(100,9) 5.37 4.56 4.33 4.25 5.3 4.01 
(100,15) 5.29 4.57 4.16 4.08 5.08 3.95 
(100,21) 5.08 4.47 4.09 3.98 4.88 3.81 
(200,8) 5.57 4.79 4.55 4.39 5.20 4.06 
(200,12) 5.43 4.66 4.43 4.35 5.13 4.03 
(200,16) 5.50 4.68 4.33 4.20 5.20 4.03 
(200,20) 5.57 4.70 4.36 4.25 5.25 4.10 


在 这 一 部 分 ， 通 过 比较 缩 边 算法 与 其 他 算法 在 相同 例子 
上 的 执行 结果 来 对 比 不 同 算法 之 间 产 生 的 解 的 效果 。 

每 个 例子 产生 规则 如 下 : n 代表 实例 中 点 的 个 数 ，R 是 
每 个 点 的 传输 半径 的 阔 值 (一 般 为 1) ,dg 是 到 处 图 Gk=(V,ER) 
的 平均 点 度 ， 可 表示 导出 图 中 边 的 稠密 度 。 对 于 每 个 给 定 的 
xS (nd), E wxw 的 方形 区 域内 随机 产生 n 个 点 ,容易 知道 
d 满足 d=nxRYw?。 对 每 个 实例 ， 比 较 每 个 算法 产生 的 导出 图 
的 干扰 度 。 随 机 产生 100 个 实例 并 用 这 100 个 实例 的 干扰 度 
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(a) n-100 时 ， 不 同 的 平均 度 4， 其 中 纵 坐 标 为 对 数 刻 度 
(a) When n = 100, the average degree D is different, in which ordinates 


are logarithmic 
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0.01 N = N - NIE ME RGA 
SE RE RE Bie 
- - - ^ - 
«EE: did: de 1 rae 
= MM 上 二 四 三 Em SCA 
aliis adem adsis qs 
0.001 SEA 
8 12 16 20 
d 


(b) 1-200 时 ， 不 同 的 平均 度 %g， 其 中 纵 坐 标 为 对 数 刻 度 


(b) When n = 200, average degree D is different, in which ordinates are 


logarithmic 

图 8 各 算法 在 不 同 随机 图 上 平均 运行 时 间 

Fig.8 Average running time of algorithms in random graphs 

从 表 2 可 以 看 出 ， 缩 边 算法 的 解 的 效果 明显 优 于 其 他 算 
法 ， 但 是 在 图 8 中 可 以 直观 地 看 出 缩 边 算法 的 运行 时 间 比 其 
他 算法 长 很 多 ， 并 且 随 着 点 对 (n,q) 的 增长 ， 运 行 时 间 呈 指数 
增长 ， 这 说 明 单 纯 的 缩 边 算法 适用 性 并 不 太 强 。 为 此 本 文 结 
合 了 增强 贪心 算法 的 策略 ， 对 缩 边 算法 进行 加 速 。 


3 — 缩 边 算法 的 加 速 


3.1 加 速 策略 
增强 型 贪心 算法 的 策略 有 很 多 ， 如 迭代 贪心 算法 已 应 用 


的 平均 值 作为 点 对 (n,q) 的 最 终 干 扰 。 根 据 下 列 的 图 8、 表 2 
可 以 直观 地 看 出 各 算法 的 解 的 效果 ， 以 及 运行 时 间 的 对 比 。 
如 表 2 所 示 ， 其 中 算 例 (x，q) 中 n 和 a 分 别 代表 相应 例子 的 
节点 数 和 平均 点 度 


T 


于 置换 流水 车 间 调 度 问 题 ，Xiong 等 人 将 随机 贪心 算法 应 用 
于 MLST 问题 0549。 这 两 类 贪心 算法 易于 实现 , 但 在 MIP 问 
题 中 随机 化 容易 破坏 图 的 连通 性 。 最 近 Carmine 等 人 提出 了 
旋转 贪心 算法 ， 易 于 实现 且 已 成 功用 于 增强 一 些 典型 问题 的 
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贪心 算法 M7。 本 文 的 第 二 章 对 旋转 贪心 算法 进行 了 解释 ， 


Ao, F: 


实验 部 分 进行 相关 了 的 实验 对 比 。 


本 文 考虑 了 结合 批量 处 理 


的 增强 策略 ， 对 缩 边 算法 进行 


了 改 ; 


机 上 批量 处 到 


TUE 


。 批 量 处 理 增强 策略 就 是 多 个 步骤 同时 运行 ， 与 计算 
基本 思想 相似 ， 这 种 策略 可 以 减少 


算法 的 运 


行 时 间 ， 使 算法 的 适用 性 更 广 。 


对 于 影响 缩 边 算法 运行 时 间 的 因素 主要 分 为 四 个 部 分 : 
a) 找 到 当前 最 大 干扰 点 ; b) 找 出 所 有 出 度 中 含有 当前 最 大 干 
扰 点 的 节点 ; c) 缩 边 并 判断 图 的 连通 性 ，q) 循 环 上 述 步 又 ， 
直至 不 存在 可 缩 边 的 节点 。 

关于 a)， 由 于 当前 连通 图 的 最 大 干扰 节点 在 经 过 一 次 缩 
边 策略 之 后 ， 会 有 全 局 性 的 改变 ， 所 以 每 次 搜索 都 需要 遍历 
全 图 。 关 于 b)， 可 以 采取 加 入 索引 的 方式 ,但 是 索引 需要 用 
到 和 矩阵 存储 ,空间 复杂 度 会 增加 ON), 并 且 在 每 次 对 节点 进 


忆 此 将 对 c) 和 d) 进 


行 缩 边 之 后 需要 对 矩阵 进行 修改 ， 难 以 对 该 部 分 进行 优化 。 
行 优化 ， 尽 可 能 减少 图 的 连通 性 的 判断 次 
数 ， 并 加 快 图 的 最 大 干扰 收敛 的 速度 。 
从 上 述 算 法 描述 中 可 以 知道 判断 图 的 连通 性 是 在 缩小 


个 节点 的 传输 半径 之 后 进行 。 通 过 大 量 实验 发 现 ， 最 大 干扰 


度 点 不 仅仅 只 


Ts 


即 需要 减少 干扰 度 的 节点 可 以 存在 多 


个 ， 但 是 在 缩 边 算法 中 上 


4 能 依次 去 对 这 些 同 样 有 具有 当前 最 大 


干扰 度 点 去 进行 缩 边 策略 ， 因 此 每 进行 一 次 缩 边 策略 就 需要 


判断 一 次 图 的 连通 性 
文 考虑 是 否 可 以 将 计算 


重 性 ， 这 样 就 会 导致 运行 时 间 过 长 。 


对 此 本 


通过 所 有 具 


再 进行 判断 
算法 2: 


几 批 量 处 理 
当前 最 大 干扰 度 点 一 次 性 的 进行 缩 边 策略 之 后 
导出 图 的 连通 性 .具体 算法 描述 如 算法 2 伪 代 码 。 


的 思想 结合 到 缩 边 算法 中 ， 


Input: 已 知 节点 集 V={1,2,.. .nm 在 平面 上 的 坐标 ， 对 于 每 个 点 


EV 的 最 大 传输 3 


Output: 图 的 传输 


连通 且 干 扰 较 小 。 


Begin 


end if 


end for 


O o N O Uu hb U NB 


10 else 


11 Xf u € Vr nax KARK 


12 end if 


do 工 maxe 当 前 导出 图 


for each 


if Vr max FH 的 点 都 
do 更 新 各 个 缩 边 点 的 


FABBRI Ri 


Ef nikns-(rnn... 


[86,1] ,满足 图 Gaz (V, Er) EE 


.rn)， 使 得 图 Gr=(V,Er) 


MSR 1€1,2...n, Vr max= 中 


vev Do 


if Iy-I max then 


do VI maxe Vr max UV 


最 大 干扰 


行 缩 边 操作 后 图 仍 连通 then 


4 
ES 


13 while I max 减 小 


14. end wh 


end 


ile 


3.2 ”加速 后 的 实验 结果 


为 了 说 明 算 法 的 改进 是 


的 缩 边 算法 


效 的 ， 本 文 对 批量 缩 边 算法 进 


行 了 对 比 实验 ， 实 验 数 据 如 表 3 和 图 9 所 示 。 其 中 算 例 (n， 


只) 中 nn 和 4d 分别 代 表 相 

根据 上 述 图 表 可 以 看 出 ， 虽 然 ] 
与 原 缩 边 算法 
优 于 其 他 算法 。 反 观 批 量 缩 边 算法 的 运行 时 间 ， 
策略 虽 不 能 将 缩 边 算法 


的 效果 对 上 


t 


应 例子 的 节点 数 和 平均 点 度 。 


比 量 缩 边 算法 的 解 的 效果 
不 大 ， 还 是 


ENS —25, (HIETH2E 
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级 到 O(02 ,但 是 其 实际 运行 时 间 却 远 远 低 于 原 缩 边 算法 的 运 
行 时 间 。 


表 3 各 算法 所 得 解 的 平均 目标 值 


Table 3 Average result of algorithms 


(nj NNA BNA RGA CGA SCA SEA  P-SEA 
(100,6) 5.15 4.56 4.39 4.36 4.80 4.08 3.89 
(100,9) 5.21 4.60 4.39 4.27 5.08 4.06 4.08 
(100,15) 5.19 4.57 4.25 4.11 5.06 3.83 4.09 
(100,21) 5.15 4.31 4.11 4.01 4.88 3.89 3.94 
(200,8) 5.58 4.94 4.76 4.80 5.22 4.11 4.20 
(200,12) — 5.53 4.68 4.38 4.36 547 4.03 4.09 
(200,16) — 5.53 4.67 4.35 4.30 5.21 3.98 4.14 
(200,20) 5.62 4.71 4.39 4.22 5:25 4.00 4.15 
1 
0.1 ENNA 
] z $ mBNA 
200 — t E pir roa 
D 9 - - - - 
- = ME = NAE 
AB NM ME 日 
MAE NAE M E Xs "CoA 
oooi Aif S REE E 
BM INEMINIBEINHERES 
- ^ - - ` - 
[RE NEM NIBE NIB 
0.0001 NY EANS = A 
6 9 15 21 


d 
(a) n=100 时 ,不 同 的 平均 度 d 其 中 纵 坐 标 为 对 数 刻 度 


(a) When n = 100, average degree D is different, in which ordinates are 


logarithmic 
1 

R z E S axNA 

0.1 E z = z 
" H E EE uL 
0.01 ` E \ E N : N E CGA 

= E SEI 日 
` z NIE NE NSESAmSCA 

NIBEM NEM VUE INE 
0.00; Wvel=d Daval NINE MRY- sea 
8 12 16 20  «P-SEA 

d 


(b) n-200 时 ,不 同 的 平均 度 d 其 中 纵 坐 标 为 对 数 刻度 


(b) When n = 200, average degree D is different, in which ordinates are 


logarithmic 
图 9 各 算法 在 不 同 随机 图 上 平均 运行 时 间 


Fig.9 Algorithms of average running time in random graphs 
+> 
4 ”结束 语 


在 以 接收 者 为 中 心 的 干扰 模型 下 ， 针 对 无 线 传感器 网 络 
中 基于 拓扑 控制 技术 的 干扰 最 小 化 问题 ， 本 文 提出 了 一 种 缩 
边 贪心 算法 。 该 算法 采用 了 已 有 算法 不 一 样 的 构图 方式 ， 通 
过 在 每 个 节点 的 传输 半径 为 最 大 的 阔 值 的 导出 图 上 采取 缩 边 
的 方式 进行 构图 在 大 多 数 随机 产生 的 算 例 上 能 产生 更 优 解 。 
针对 其 运行 时 间 过 长 这 一 缺陷 ， 本 文 将 批量 处 
缩 边 算法 中 ， 大 大 缩短 了 其 运行 时 间 ， 相 较 于 其 他 算法 使 得 


批量 的 加 速 


的 时 间 复 杂 度 从 CCD 降低 一 个 数量 


其 能 在 几乎 相同 的 时 间 内 获得 更 优 的 结果 。 


定稿 LIE RE 
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